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Based upon mesoscale simulations of binary mixtures with very low surface tension and positive 
disjoining pressure (frustration), we measure the correlation length of the stress field within the 
flowing mixture, as a function of the frequency of the applied load. Two scaling regimes are clearly 
identified, with a sharp crossover between the liquid and solid regions, with exponent —1/2 and —I, 
respectively. The crossover correlation length is shown to be closely connected to the typical scale 
of the coherent excitations of the density field, whose collective dynamics is shown to be responsible 
for non-trivial rheological effects including the emergence of yield stress. 
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INTRODUCTION 



Soft-Glassy (SG) materials exhibit features of both liq- 
uid and solid states of matter [l], and consequently 
they present a major challenge to standard concepts and 
tools of thermod yna mics and non equilibrium statisti- 
cal mechanics p- 18|. Despite their disparate physico- 
chemical nature, such materials share several generic fea- 
tures, primarily an anomalous relaxation time to equilib- 
rium and a non-linear relation between the applied stress 
and the resulting strain rate (non- Newtonian rheology). 
This behaviour is sometimes accompanied by a so called 
yield-stress, i.e. a minimum stress threshold, below which 
the material does not flow, showing solid-like behaviour 
Q • The microscopic roots of such an intriguing rheology 
have made the object of intense scrutiny in the recent 
years, via molecular dynamics simulations [17|, meso- 
scopic models 0, E3, 0, as wen as experimental work 
0, S H 16 1 • Although much has been learned thanks 



to these efforts, a commonly agreed theoretical picture is 
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yet to emerge 

An important question which underlies the aforemen- 
tioned phenomenology concerns the basic mechanisms 
controlling the crossover from liquid to solid-like be- 
haviour of SG materials. In particular, in this Letter, we 
exhibit evidence that there exists a dynamic crossover 
scale, Idcs, marking the transition between the two 
regimes, i.e., phenomena correlating below (above) Idcs 
are classified as liquids (solids). We propose a well- 
defined operational procedure to measure this character- 
istic length-scale, based on the spatial correlation of the 
velocity (stress) field. It is reasonable to expect that such 
correlation length-scale bears a link with other length- 
scales discussed in the literature, namely the shear trans- 
formation zones [l9| and cooperative length-scales Q . In 
our work we employ a mesoscopic model of binary fluid, 
where a suitable frustration mechanism at the interface 
is introduced by combining a small positive surface ten- 
sion, promoting highly complex interfaces, with a posi- 
tive disjoining pressure, inhibiting interface coalescence. 



The mesoscopic kinetic model considers a binary mixture 
of fluids A and B, each described by a discrete kinetic 
Boltzmann distribution function f s (r,v~l,t), measuring 
the probability of finding a representative particle of fluid 
s = A, B at position r and time t, with discrete velocity 
Vi, where the index i runs over the nearest and next-to- 
nearest neighbors of f in a regular two-dimensional lat- 
tice [2(| • By definition, the mesoscale particle represents 
all molecules contained in a unit cell of the lattice. The 
distribution functions of the two fluids evolve under the 
effect of free-streaming and local two-body collisions, de- 
scribed by a momentum-relaxation to a local equilibrium 
(/ ) on a time scale r for both fluids (s = A, B): 

f s (f+ yl,vl, t + 1)- f s (r, v l7 1) = 

- \ (/. - /i e9) ) (?, vl;t) + F s (r, v-;t). 

The term F s (r, vl, t) includes a variety of interparticle 
forces. First, a repulsive force with strength parame- 
ter Qab between the two fluids ensures phase separation 
[20j (see triangles in figure [1]). Furthermore, both flu- 
ids are also subject to short range self-attraction (con- 
trolled by strength parameters Qaa,i < 0, Qbb.i < 0), 
and mid-range self-repulsion (regulated by strength pa- 
rameters Qaa,2 > 0, Qbb,2 > 0). In this way, one has the 
basic ingredients for ideal phase separation plus forces 
within each species whose role is twofold. First, they set 
the desired non-ideal (excess) pressure in the bulk. Sec- 
ond, they introduce an effective "long-range" repulsion 
between the two species which frustrates phase separa- 
tion. Despite their inherent microscopic simplicity, such 
dynamic rules are able to promote a host of non-trivial 
collective effects [2(| Hi]]. Most notably, whenever the 
self-specie attraction (Q ss ,i) and repulsion effects (Q ss ,2) 
come close to an exact balance, two distinctive phenom- 
ena take place: low surface tension is achieved, thus 
producing a large interface between the two fluids, and 
a positive disjoining pressure, inhibiting coalescence of 
the bulk phases. This is what we call frustration in our 
model. The emergence of a positive disjoining pressure 



FIG. 1: The order parameter of a thin film composed of two 
non ideal interfaces separated by a distance h. Each interface 
is characterized by a surface tension T. The order parameter 
is defined as (f> = ±1 in the A(B) bulk phase and = 
at the interface location. Two cases are considered: (a) a 
mixture with purely repulsive fluids; (b) the same mixture 
plus a frustration mechanism. The presence of frustration 
at the interface translates into a positive disjoining pressure 
11(h) as the film gets thinner (see inset). 



II(/i) is shown in figure [IJ where we consider a thin film 
with two non ideal interfaces, separated by the distance 
h [28| . Following Bergeron [22[ , wc write the relation for 
the corresponding tensions following the Gibbs-Duhem 
relation: 



f n(h) 

T f (h) = 2T+ h 

Jn(h=oo) 



dn 



(1) 



where T is the bulk value of the surface tension and Tf 
is the overall film tension, whose expression is known 
in terms of the mismatch between the normal and tan- 



gential components of the pressure tensor [23|, |24| , Tf = 
f_°°(Pjf — Pr{x))dx, where, in our model, Pn~Pt(x) = 
p s (x) can be computed analytically [IBf . With the knowl- 
edge of the relation s(h) = Tf(h) — 2T, the computation 
of the disjoining pressure is then straightforward: a sim- 
ple differentiation of s(h) allows us to compute the first 
derivative of the disjoining pressure, ds jl^ = /i4r- This 



dh 



information, supplemented with the boundary condition 
U(h — > oo) = 0, allows us to completely determine the 
disjoining pressure of the film (see figure [I}. The diffuse 
interface profile emerges spontaneously from the underly- 
ing mcsoscopic dynamics, once the strength of the under- 
lying interactions is properly tuned. We wish to remark 
that the humped-structure at the frustrated interface (see 
figure [IJ would disappear by lifting the tight competition 
between short-range attraction and mid-range repulsion, 
thus showing that frustration is indeed essential to re- 
produce the correct interface physics. In the absence 
of such a stabilization mechanism, the rheology of the 
binary mixture would reduce to a plain Newtonian be- 
haviour. 

In the following, we turn to the analysis of the basic 
rhcological behaviour of two dimensional emulsion and 



FIG. 2: Typical density configurations in our numerical simu- 
lations. We simulate cases with 50:50 (case (a), upper panel) 
and 30:70 (case (b), lower panel) mass ratio between fluids A 
and B, respectively. Case (a): we simulate islands of fluid A 
embedded in fluid B, all entrapped within a "vesicle" confined 
by a belt of fluid A. The initial condition is a periodic per- 
turbation of the uniform background density, with wavelength 
lo — 64 and zero flow. Case (b): we prepare the system with 
a collection of polydisperse "packings" of fluid A (B), sepa- 
rated by layers of fluid B and being stabilized by a positive 
disjoining pressure (see also figure [TJ • In both cases, the A 
rich (B rich) region has approximately pa = 1-2, ps — 0.2 
(pA = 0.2, ps = 1.2). The hydrodynamic viscous ratio be- 
tween A rich and B rich regions is 1. The coupling parameters 
are Qaa,\ = -9.0, Qaa,i = 8.1 Qbb,i = -8.0, Qbb,i = 7.1, 
Gab = 0.587 and p = 0.83 [IJ. 



foam-like fluids, such as those described in figure We 
perform two different experiments: in the first case (the 
"full periodic" experiment) the fluid is placed in a square 
box L x x L y = 512x512 with full periodic boundary con- 
ditions. A forcing is then applied in the stream-flow di- 
rection, say x, with the property F x {y) = -Fbsin f^p^Jj 
where the constant F$ is tuned so as to reproduce for 
a stationary fully Newtonian case the velocity profile 

u x{y) = ?7osin {^j^j- When the frustration mechanism 
is switched on, the velocity profile deviates from the New- 
tonian prediction. From the averaged velocity in time 
and in the stream-flow direction we are able to extract 
the local strain rate, while the stress is provided by the 
integral of F x {y)- Such local rheological measurement is 
displayed in figure |3] (circles) and clearly indicates the 
existence of a yield stress try ~ 10~ 4 . This result is con- 
sistent with the rheology obtained in a second experiment 
(the "Couette" experiment) where the system is confined 
between two walls (y = 0, y = L y ) with periodic inlet- 
outlet boundary conditions. The upper boundary is set 
to a constant velocity, Uw, resulting in a constant aver- 
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aged stress. The stress-strain rate plot given in figure [3] 
(squares) reveals that both rheological measurements are 
in a fairly good agreement. We now turn our attention to 
the case of a sinusoidal strain applied in the confined Cou- 
ette experiment. We increase the peak strain amplitude 7 
and compute the peak stress ap, while keeping constant 
the oscillation frequency uj = 1 CP 4 , in line with the anal- 
ysis reported by Mason & coworkers Q. The results are 
shown in the upper inset of figure [3] A power law in the 
peak stress-strain relation, ap ~ 7, is observed for the 
lower strains. For large 7, the system starts to yield and 
the peak stress exhibits a sub- linear power law behaviour, 
as observed in 0] and the yield stress ay is, within a fac- 
tor 2, consistent with the one obtained in the main plot of 
figure [3l Next, we investigate the response of the system 
when the upper boundary is oscillating according to the 
relation uw(t) = Uw cos(ujt), with uj the load frequency 
and Uw = 0.001. Notice that the strain is now function 
of uj, i.e. 7 ~ XJ-wju). In this setup, the stress is typi- 
cally expressed as a(t) = G'{uj) sin(u;i) + G"(oj) cos(wt), 
where G' is the storage modulus and G" is the loss mod- 
ulus. In the bottom insert of figure [3] we show G' and 
G" as a function of uj. In typical SG materials, say con- 
centrated emulsions, there is a prevalence of G", up to a 
crossover frequency, where the fluid component G" takes 
over. It is apparent, looking at the inset, that G'(uj) 
(solid response) prevails over G"{uj) (liquid response) at 
all frequencies up to w ~ 5 - 8 x 10~ 3 . 
The main goal of the present letter is to investigate the 
physics of spatial structures responsible for the complex 
and non linear rheology described in figure [3] To this 
purpose we consider the spatial stress correlator: 



Ca(y) 



(a{L y ,t) a{y,t)) 
(o*{L Vi t))V*(o*{v,t))W 



(2) 



where a(y,t) denotes the averaged stress in the stream- 
flow direction and the brackets denote time average (typ- 
ically over 10 6 time-steps). In a liquid system with vis- 
cosity i/jv, the fluid velocity u x {y,t), averaged in the 
stream-flow direction, oscillates in time with frequency 
oj and decays as e~( Lv ~ v '' lf , where If ~ \J i>n/u. As 
a consequence, the stress a{y, t) shows the same behav- 
ior and the stress correlation ((2|) behaves as C a (y) = 
cos(2n(L y — y)/lf). In a solid-elastic system, we expect 
C a {.y) = cos(2n(L y - y)/l s ), where l s cx c s /uj, c s being 
the sound speed. Thus we conclude that, in our system, 
as well as in many other SG systems, the correlator can 
be written as C„(y) = cos(27r(ij / — y)/l{uj)). As a result, 
we can use the scaling 1(uj) vs uj to quantitatively assess 
the transition from liquid- to solid-like behavior. Notice 
that the computation of 1(uj) can be carried out only for 
large enough L yi i.e. when inertial effects become im- 
portant. In figure IH we report the results for 1(oj), as 
extracted from the spatial stress correlator C a (y) in a 
domain of size L x = 1024 and L y = 1024 for the setup 
of case (a) of figure [5] (results related to the setup of 



case (b) are presented in figure [5]). The wall speed is 
fixed to Uw — 0.001, whereas the load frequency is var- 
ied in the range 10 -4 < oj < 10~ 2 . The simulation data 
provide a very neat indication of the two regimes, with 
exponents —1/2 and —1, respectively. Moreover, a sharp 
cross-over between the two regimes is clearly observed at 
a crossover frequency uj c w 3 x 10 -3 , which identifies 
a correlation length, Idcs ~ 100, in lattice units. This 
length is of the order of a few l (see caption of figure 
[2]) , i.e. a few times the characteristic size of the density 
wavelength in the initial condition, indicating that the 
collective excitations of the density field in the solid-like 
phase are basically lumps of fluid within a given belt, 
surrounded by their nearest neighbors. We also repeated 
the numerical simulations for the case of a fluid with- 
out liquid-liquid repulsion Gab = 0, thus recovering the 
Newtonian behaviour for a homogeneous fluid. Consis- 
tently with expectations, we note that the liquid scaling 
(squares) extends all the way into the low frequency re- 
gion. A possible way to rephrase the results of figure 
IH is to define a frequency-dependent effective viscosity, 
v e ff{uj) = oj1 2 {oj) « c 2 /uj, that highlights the divergence 
of the viscosity in the low- frequency limit uj — > 0. In 
the inset of figure |4] we show v e f f(oj) / Vjv as a function of 
oj, where corresponds to the Newtonian viscosity ob- 
tained for Qab = 0. In the same inset, we plot the peak 
stress ap normalized with respect to the corresponding 
Newtonian value. Both v e ff and ap increase as 1/oj for 
small uj, which is consistent with the interpretation of 
uj1 2 {uj) as an effective viscosity. Let us remark that, in 
a Newtonian liquid, ap ~ vjyUw / Ly should be indepen- 
dent of uj whereas in a solid-clastic system ap is propor- 
tional to 1/uj. The picture emerging from our results is 
that the fluid behaves like a liquid at short-scales, where 
the leading equilibration mechanism is the local relax- 
ation in momentum space, and like an amorphous solid at 
large scales, where the correlated dynamics of the collec- 
tive- degrees of freedom undergoes anomalous slow-down, 
with an ensuing substantial increase of the viscosity (slow 
relaxation in configuration space). 

An important point is to investigate whether the no- 
tion of Idcs may be related to the existence of yield stress 
cry. Following [3j and our interpretation of Idcs, we ar- 
gue that ay should be given by the ratio of the surface 
tension V and Idcs- in our case, T ~ 10~ 2 suggesting 
ay ~ 10 -4 in agreement with the basic rheological results 
shown in figure |3] To gain a deeper physical insight, we 
consider a bi-liquid emulsion system with monodisperse 
distribution of droplets sizes (see case (b) in figure [2]) 
within a confined Couette geometry with upper wall ve- 
locity Uw constant in time. Next, we perform different 
experiments by varying the nominal shear S = Uw/L y 
and compute the average stress (a) and the stress corre- 
lation function ([2]) at the channel center y = L y /2. In the 
inset of figure [SJ we show three cases of C a (y) obtained 
for different values of S. For each C a (y) we can estimate 
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FIG. 3: Basic rheological properties for the soft-glassy con- 
figuration of case (a) in Figure (2J see text for details. In the 
upper inset we report the peak stress op produced by ap- 
plying a sinusoidal strain with peak amplitude 7 at a given 
frequency ui = 1CF 4 . Shown in the lower inset, is the be- 
haviour of G' {uj) and G"(u), where G' is the storage modulus 
and G" is the loss modulus. 
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FIG. 4: Stress correlation length as a function of the load 
frequency. The two straight lines refer to slopes —1 (solid 
behavior) and —1/2 (liquid behavior). The dynamic cross 
over length is estimated to be Idcs ~ 100. The squares refer 
to the Newtonian case (i.e. Qab = 0) which follows the slope 
— 1/2. Inset: we report the effective viscosity u e ff = uil 2 (u) 
and the peak stress op normalized to their Newtonian values. 
For small ui, both quantities show the behavior l/u. 



the decay length Id as C a (L y /2 ± Id) 



Finally, in 



the main plot of figure we show (a) as a function of 
Id 1 1 DCS- If our interpretation of ay ~ T /Idcs is cor- 
rect, we expect that (a) jay > 1 for Id ~ Idcs, which is 
exactly what we observe. The physical meaning of our re- 
sults implies that the condition required to set the system 
in motion (i.e. for (er) > cry) is that the characteristic 
fluctuations should be correlated at relatively small scale, 
i.e. smaller than Idcs- I n summary, the physical mean- 
ing oil dcs appears to be consistent with the solid/liquid 
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FIG. 5: Average stress (normalized with the yield stress) ver- 
sus the decay length Id normalized with Idcs- Data refer to 
a Couette experiment (see also figure [3} at different nomi- 
nal shears S = Uw/Ly for a bi- liquid emulsion system with 
monodisperse distribution of droplets sizes (see case (b) in 
figure [5]). The inset shows the profiles of the stress correla- 
tor for three cases: Uw = 0.005 (circles), Uw = 0.01 (stars), 
Uw = 0.04 (squares), while L y is kept fixed. 



transition and the existence of a yield stress. Based on 
the observations prompted by mcsoscale simulations of 
a binary mixture with competing attractive/repulsive in- 
teractions, in this paper we have proposed an operational 
procedure to identify a typical scale separating fluid-like 
and solid-like behaviour. Our approach can be tested 
in laboratory experiments, similar to those recently dis- 
cussed in [121 ] . Our model provides a remarkable blend 
of conceptual simplicity (all interactions are prescribed 
in terms of elementary lattice interactions) and compu- 
tational efficiency (2t| . 
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